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SUMMARY 

The usual approxLmations made in calculating the pressure 
changes of a compresaihle fluid floving through heat -exchanger 
passages introduce appreciable errors in the range of hig^-flow 
Mach numbers and high rates of heating. Existing methods for 
obtaining accurate results over the entire Mach number and rate- 
of-heating range require numerical integration for each specific 
set of conditions and are therefore too tedious for general 
application. 

In the present paper ein analysis is made of the compreesible- 
flow variations occurring in heat -exchanger passages. The results 
of the analysis describe the flow and heating characteristics for 
which specific flow passages can be treated as segments of gen- 
eralized flow systems. The graphical representation of the flow 
variations in the generalized flow systems can then be utilized 
as working charts to detennine directly the pressure changes 
occurring in any specific flow passage. On the basis of those 
results, vorkliag charts are constructed to handle the case of air 
heated an constant wall temperature under turbulent-flow condi- 
tions. A method is given of incorporating the effect on the heat- 
exchanger flow process of high temperature differential between 
passage v 11 and fluid as based on recent NACA experimental data. 
Good agreement is obtained between the experimental and the chart 
pressure-drop values for passage-wall average temperatures as high 
as 1752° R (experimental limit) and for flow Mach numbers ranging 
from 0.32 to 1.00 (choke) at the passage exit. 
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IBTBODDCTION 

Tho rational design, and calculation of performance of aircraft 
heat exchangers wherein heat is added to or subtracted from a high- 
speed compressible fluid stream requires not only the knowledge of 
the appropriate heat-transfer and friction-drag coefficients but 
also the mathematical means of utilizing these coefficients to 
describe accxorately the pressure and temperature variations along 
the length of the heat-exchanger flow passage. 

The basic differential equation describing the pressure varia- 
tions of a compressible fluid under the simultaneous action of 
friction and heating or cooling is well known for the ideal case 
of one-dimensional f lowj however, no exact closed-foim solution of 
this differential-flow equation has yet been obtained even for 
specially chosen boundary conditions. Simplified solutions have 
been evolved (references 1 and 2) by resorting to approximations 
that are sufficiently valid for flow at relatively low speeds 
(corresponding to Mach nimibers of the order of 0.4 and less) and 
for moderate rates of heat input to the fluid (corresponding to 
temperature differentials between wall and fluid of the order of 
300° R and less). However, as Illustrated by the experimental 
results of reference 3, the eirors introduced by the approximations 
in calculating the pressure variations along a flow passage 
Increase so rapidly with increase of Mach numiber and rate of heat 
input tha t for many heat -exchanger problems the simplified solu- 
tions are, at the most, only rou^ approximations. 

Although the differential-flow equation is not susceptible to 
formal Integration, it is readily amenable to numerical methods. 
Methods for numerically integrating the differential-flow equation 
for specific heat- exchanger conditions are discussed in refer- 
ences 4 ft rid 5. Although the methods devised reduce the labor 
involved in obtaining numerical solutions it la desirable in the 
calculation of pressure drop to obviate the necessity for per- 
fozming a numerical integration for each set of conditions of heat- 
exchanger operation. 

In the present paper, an analysis of the heat -exchanger flow 
process is made that describes the generalization conditions for 
which specific flow passages can be treated as segments of gen- 
eralized flow systems. The results of the analysis provide the 
basis for the construction of working charts that enable deter- 
mination, without individual Integration, of the pressure variations 
of a compressible fluid flowing through heat- exchanger passages. 
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Working charts are herein presented for the specific case of 
air floving in tiarhulent motion throxigh smooth constant-area passages 
wherein heat is added to the air stream hy the passage walls which 
are at a constfint temperature throughout their length. A constant 
ratio of the specific heats for air eqvial to 1.400 is used in the 
preparation of these charts. The range of variables covered in the 
charts are: (a) ratios of passage-wall temperature to air temper- 

ature (absolute) from 5.00 to 1.11, (b) total-momentum parameters 
equivalent to entrance Mach numbers from 0.15 to 0.80, and (c) pas- 
sage length-diameter ratios from zero to a value sufficient in most 
cases to produce choking at the jassage exit. 

Althoxi^ in the construction of the charts the heat-transfer 
coefficient and friction factor are expressed by the standard 
turbulent-flow equations used in present-day heat -exchanger work, 
a method is given for incorporating the effects of hl^ temperature 
differential between passage wall and flxild on heat-transfer coef- 
ficient aJid friction factor as indicated by the recent experimental 
data of reference 6. A method of incorporating the e:^erlmental 
res\d.ts of reference 6 in the heat -transfer ralatlon used for 
temperature -rise calculation is also given. Calculated and measured 
values of pressure drop and tenperature rise are conpared. 

An example is presented that illustrates in detail the method 
of using the working charts and of using the results of reference 6 
in solution of a typical problem involving flow through a heat- 
exchanger passage. Application of the method of reference 5 for 
calculating the pressure changes of a compressible fluid flowing 
in a passage with an abrapt Increase of cross section (as obtained 
at exit of heat exchanger) is discussed and illustrated. 

Although of a preliminary nature, this paper is presented to 
fulfill the current urgent needs for this type of Infoimation. 


SyMBOLS 


A cross-sectional area of flow passage, (sq ft) 

Cp specific heat of fluid at constant pressure, (Btu/(lb)(°E)) 
Dq equivalent diameter of flow passage. 

Dp drag force due to friction, (lb) 
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friction factor. 



free -flow airea ratio 

mass conversion factor, 32.2, (ib/slvig) 

heat-transfer coefficient "between vail and fluid, 
(Btu/(sec)(sq ft)(9E)) 


ratio of actual value of friction factor to that given "by the 
standard expression, 0.046 

ratio of actual value of heat-transfer coefficient to that 

KDq 

given "by the standard expression, = 0,023 

thermal conductivity of fluid, (Btu/(sec)(sq ft)(°E)/(ft)) 
passage flow length, (ft) 

Mach number 

mass flow of fluid, (slugs/sec) 
total pressure, (l"b/sq ft a"bsolute) 
static pressure, (Ih/sq ft a"bsolute) 

gas constant for fluid, (for air, equal to 53.35), (ft-lb/(l'b) (°E) ) 

wetted perimeter, (ft) 

total telirperature of fluid, (9 r) 

teB 5 »erature of passage walls, (^) 

static ten?)erature of fluid, (9 r) 

fluid velocity in flow passage, (ft) /(sec) 
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X downstream distance from reference station in flow system, (ft) 

Z parameter determined from tatle on charts 

y i^tio of specific heats at total temperature of fluid (asstuned 
constant as 1.400 in air charts) 

(i viscosity of fluid, (lh)/(ft)(sec) 

Up viscosity of air at 519° R, 12.3 X 10"®, (lh)/(ft)(sec) 

p mass density of fluid, ( slugs /cu ft) 

Subscripts : 

av flviid conditions evaluated at average fluid temperature in flow 
passage 

en entrance of given flow passage 

ex exit of given flow passage 

w fluid conditions evaluated at average wall ten 5 >erature of flow 
passage 

X any station in flow system 

0 reference station, defined as station in flow system at which 
ratio is equal to any of values 5.0, 3.0, 2.0, or 1.5 

Die following groupings of variables are involved: 


aiV+pA 

total momentum 

aST+pA 

total-momentum parameter 

m/jgE?T 


jpA 

static -pressure parameter 

m YgElT 


PA 

total -pressure parameter 

m-^ gE?T 
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V 


velocity parameter 




0.2 


effective length- diameter 
ratio as measured from 
reference station to any 
station in flow 
system 



effective length- diameter 
ratio of flow passcige 


AMLYSIS 


The steady-flow process occurring within a constant-area heat- 
exchanger passage is one involving the simultaneous action of fluid 
friction and heat transfer. A general analysis is herein made of 
the compressible-flow variations obtained during this flow process. 

It is shown that, by appropriate limitations of the flow arifl 
heating conditions, individual flow processes become special cases 
of generalized flow systems; as a result generalized graiiiical 
representation of flow systems can be constructed that are applicable 
for the direct deteimination of the flow variations in a large 
number of different flow passages. The results of the analysis are 
then utilized for the construction of charts for air heated at 
constant wall temperature, assuming the standard heat- transfer and 
friction-factor eq^\iatlons for turbulent flow through smooth pas- 
sages. A method of incorporating the effects of high temperature 
differential between passage wall and fluid on the flow process, as 
indicated by the recent experimental data of reference 6, is pre- 
sented and checked with experiment. 


For an over-all evaluation of heat -exchanger perfommice, the 
flow variations at the heat -exchanger entrance due to the flow-area 
contraction and at the heat-exchanger exit due to the flow-area 
enlargement must also be considered. Inasmuch as the entrance 
losses aire usually small compared with the over-all losses across 
the heat exchanger, the flow process at the hsat-eicbenger entrance 
can generally be assumed isentropic. For a sharp-edge entrance, 
a more accurate determination of the pressure change in the flow 
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proceaa at the heat-exchaiiger entrance can be obtained by aubtractlng 
the presaure drop due to the vena-contracta, aa based on Incompreaslble- 
flow considerations, from the pressure calculated for compressible 
Isentroplc flow. The losses at the heat-exchanger exit are usually 
of appreciable magnitude and requlire accurate eTaluatlon. For this 
purpose, the analysis presented In reference 5 of the compressible- 
flow variations across a sudden enlargement Is herein applied for the 
case of subsonic flow that Is generally encountered In heat- exchanger 
practice . 

Description of flow process. - One foim of the dlfferentlal- 
momentxam equation describing the one-dlmenslonal steady- state motion 
of a compressible fluid in a constemt-area passage under the combined 
Influence of friction and heat transfer Is: 

dCmV+pA) + dDj> = 0 (l) 

Prom the conservation of energy and mass equations and the perfect - 
gas law, reference 5 shows that the total-momentum parameter 

i 0 uniquely related to each of the flow parameters M, — — — , 

mV^ 

— ^ , and ^ for any value of y associated with the total tem- 

P * 

perature of the fluid. Knowledge of the variations of and 

m/^gRr 

T during the flow process Is therefore sufficient to completely 
specify the variations of all the fluid-flow conditions (p, P, V, 
and t) . 

The variation of the total-momentum parameter Is due 

m VgKT 

to the variations of both the total-momentum mV+pA and the total 
temperature T. Differentiation of the total-momentum parameter 
with respect to these two variables for constant mass flow and area 
gives : 


i ( _ d.(my+pA) _ 1 mV+pA dT 

\m-^gBT/ m VgET ^ m VfiiOT ^ 

The differential drag force dDp Is 

dDp. = F i pV^A ^ 

2 Pq 


( 2 ) 


( 3 ) 
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or, more conveniently, 


eiDp 


2FMVa(i) 


(4) 


The variation of total temperature T 1 b given hy the differential 
relation equating the heat transferred from the wall to the fluid 
to the heat absorbed by the fluid. 

mgCp dT = hs (T^-T) dx (5) 


The temperature recovery factor in the boundary layer is assvnned to 
be unity in equation (5). Negligible errors are introduced by this 
assumption for gaseous fluids and for the practical range of tem- 
peratiire differentials used in beat exchangers. 


Prandtl's extension of Reynold’s analogy between fluid 
friction eind heat transfer for flow through smooth-wall passages 


gives 


F 

2 


as a function principally of the Prandtl number 




In view of this relation, equation (5) is rearranged and presented 
as 


dT 



(T^T) F 4 (i) 


(6) 


From equations (1), (2), (4), and (6), the variation of the total- 
momentum parameter during the flow process is expressible as; 


^ /mV4£^\ 

VmV^y 


dT 


1 mV+pA dT 


/_ 

/ C 


m 


pPgV 

^ j 

2 / 




( 7 ) 


(Tw-T) 


As previously pointed out, 


m7+pA 

m.^gRT 


in the maimer dictated by the laws 
mass, the perfect-gas law, a n d the 


V 

is related to ^ 

of conservation of energy and 
relation between y and T; 
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h 

equation (7) thiie InTolyee only the Tariablee , T^, 

m ^gKr ^ 

and T. The conclusion can then he dravn that for a giyen fluid the 


yariation of 


aV+pA 

a VeiRT 


during the flow process is uniquely deteiained 


hy the yariation of T duidng the process pTOyided that 


CpPgV 

F 


and 


2 

Ty are constant or are expressible as functions of only mV +pA 

V tiA “ 

and T. Because ^ ta and T coapletely define all the fluid- 
B VgBT 


flow conditions, the parallel generalization can be aade that the 
flow yarlatlons of a fluid along a heat- exchanger passage are 
uniquely detemlned by the yariation of T along the passage when 


the local yalues of 


h 

and T|f are constant or are expressible 

2 


CpPgY 

F 


as functions of only the local fluid-flow conditions. In partic- 
ular, the flow yarlatlons are then explicitly independent of the 
position yarlable x. Thus when the generalization conditions are 
satisfied, any station in a giyen flow system can be considered as 
the initial starting point of a real flow process; conyersely, 
Indlyldual flow processes can be considered as special cases of 
generalized flow systems. 


The results of integration of equations (6) and (7) for a num- 
ber of generalized flow systems can be presented in chart form for 
the direct detezalnation of the flow yarlatlons in a large number 
of indlyldual flow passages when the generalization conditions are 
satisfied. This fact is used as the basis for the construction of 
the charts presented herein. 


P3ressure-drop charts for air heated at constant wall tempera- 
ture. - Working charts were constxoicted on the basis of the fore- 
going analysis to enable deteimination, without indlyldual integra- 
tion, of the pressure yarlatlons sustained by air flowing through 
constant-area passages irtierein heat is added to the air at constant 
passage-wall temperature. Conyentional heat-transfer and friction 
relations were used in the preparation of the charts. Subsequent 
to the completion of the charts, additional data were obtained 
(reference 6) that Indicate an effect of passage-wall temperature 
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on heat-transfer coedET iclent and friction factor that hecomes 
Important at hl^ temperature dldfferentlals hetveen vail anH fluid. 
Hence, and Ep are herein Introduced as correction factors 

that account for departure of the actual heat-transfer and friction 
phenomena from the phenomena expressed hy the conventional relations. 
The method of evaluating and ^ in accordance vlth the recent 
data of reference 6 is reserved for a later section. 

The relations and assumptions used In the Integration of equa- 
tions (6) and (7) are as follows: 

(a) The ratio of the specific heats y Is taken constant 
as 1.400 

(b) The relation for friction factor Is taken as 


When Ep = 1.00, equation (8) reduces to the conventional relation 
(reference 7, p. 119, equation (9a)). Inaeonuch as (x enters In F 
only to the 0.2 power, It Is assumed constant In the Integration of 
equation (6). Large variations In |x are accounted for In the use 
of the results by taking an average value of n in the flow pas- 
sage (that is, IX = wt£iy) . The factor % Is also assumed constant 
in the Integration of equation (6) aM, as indicated later, is 
evaluated for the average fluid conditions existing in the flow 
passage. 

(c) The relation for heat-transfer coefficient Is taken as 


When Eh = 1,00, equation (9) reduces to the conventional relation 
(reference 7, p. 168, equation (4c)). 

From equations (8) and (9), 



( 8 ) 



( 9 ) 



( 10 ) 
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The simplifications are made that: first, (^) equals 1.186 

which represents, within 1.5 percent, the value for air over the 
extreme temperature range 460° to 2500° E; and second, ~ is 

Kji 

equal to unity in accordance with Prandtl's extension of Reynold's 


h 


I rm. 

analogy. Thus, — — is taken constant as 1.186 in the inte- 

2 

gratlon of equations (6) and (7). 


Two modes of heat addition of present interest are character- 
ized by having (a) constant passage-wall temperature (T^ = a con- 
stant) and (b) constant rate of heat input along the passage 
Isngth |[h(Ty-T) = a constant J . The generalization conditions are 
satisfied for both these modes of heat addition for the foregoing 
assumed heat-transfer and friction phenomena during the process; 
generalization is therefore possible. Only the case of constant 
passage-wall temperature is specifically treated herein. 


For convenience of analysis and of presentation of results, 
the differential variable in equations (6) and (7) is changed from 

T to In the integration of equations (6) and (7), the lower 

integration limits are designated by the subscript 0. When equa- 
tion (8) is substituted for F in equation (6) and p is considered 
to be equal to |iav^ "the parameter 



'eff 



(pVDe)"°-^ Kp ~ 


( 11 ) 


is introduced. This parameter, which is proportional to the product 
of the actual friction factor FKp. and the length-diameter ratio 

is of the same order of magnitude as 
i'e ^e 


On the basis of the foregoing manipulations and assumptions, 
equation (6) can be formally integrated to give: 



( 12 ) 
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where 


^ = 1 + 


\ T 


- 1 


1 - e 


•0.00568 


'e/eff 


(13) 


Integration of equation (7) is nronerically performed to give 
m7+pA 

m <\/gET 

for each assumed initial pair of values of 


values of 


in one to one correspondence with values of 

^ mV+pA ^ 


T /o' 

Inasmuch as the flow variations are explicitly independent of the 


m VgRT/o 


and [ 


w 

T 


position variable 
Bnd 


T„ 

each corresponding pair of values of — 


__ obtained during a single integration integration for 

“VeBT \ /T„.n ~ L 

and 


one pair of values of f . 

\m V ^/0 


T /o 


can be considered as 


the entrance conditions of a flow passage; thus each integration 
gives results for an infinite number of flow passages whose entrance 

conditions are restricted to the one-to-one relation between 

m 

and obtained in the integration. All possible combinations of 
1 T 

mY+|A w ^ encountered in heat -exchanger practice 

m VgKP 1 ^ “o i:- 

are obtained if the integration of equation (7) is carried out for a 
single value of ( - 7 ^ ) and a range of values of f , How- 

^ T /Q \m Jo 

ever, in order to obtain Increased accuracy for the entire range of 

Interest, the integration was actually carried out and is presented 

'T„\ 

From the results of numerical Integra- 


for four values of 


3 

T 


tion of equation (7) and from equations (12) and (13), the values of 


mY+pA 


m VgBT 
each value of 


in one-to-one correspondence with values of ( ^ 


and 


obtained. 

later. 


/ mY+pA 


for 


e /ef f 


used in the integration are 


The plotting of these insults in chart form is described 


It is convenient in the use of the charts to recast equa- 
tions (12) and (13) as follows: 

Solution for (• 7^1 from equation (13) and substitution for 
T V^e/eff 

7 =p from equation ( 12 ) gives for x = in which case sub- 

^0 

script X is replaced by subscript en, 
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(14) 


Eq.uatlon (13) Is rewritten by substituting subscripts ex 
en for x and 0, respectively, in wblch case x = L 


Evaluation of and Kp. - The preliminary heat- transfer 

and pressure-drop data of reference 6 are used to determine the 
average values of Kjj and % that account for the effect of high 
temperature differentials between passage wall and fluid on the 
heat-transfer and pressure-loss phenomena in smooth-flow passages. 
These data were recently obtained at the NACA Cleveland laboratory 
in an Investigation conducted with air flowing through a tube that 
was electrically heated to average wall temperatures of from 710° 
to 1700° E. Eeference 6 shows that conflation of the heat-treuisfer 
data according to equation (9) wherein was taken as unity and 

the physical properties Cp, jji, p, and k of the fluid were 

evaluated at the average bulk temperature if suited in a separation 
of the data as the temperature level of the wall varied. It was 
foiind, however, that satisfactory correlation of the heat-transfer 
data over the entire range of wall temperature is obtained by the 
equation: 



1 1 - e 


-0.00568 


(^eXff (15) 



(16) 


\dierein 



( 17 ) 


Equation (16) can be rearranged as follows: 
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If 1 b asBinned constant for air and if the slurplification Is made 

that p is inversely proportional to T rather than t, the 
average value of Bh in a flow passage is obtained from equa- 
tions ( 9 ) and (18) as: 


Kh 


.022 /^^av\ * 
0.023 \T^J 


*P,ay 



(19) 


The average values of T, li, and Cp used in equation (19) 
are consistent with the deteimination of an average value of 
for the flow passage. 


Die pressure -drop data of reference 6 are presently in the 
analysis stage; no definite conclusions have been reached regarding 
the effect of hl^ tenq)erature differentials between passage wall 
and fluid on friction factor. 3ii the absence of such information 

h 


it can be tentatively assumed that 


CpPgV 

— remains substantially 


2 

consteint with variation of wall temperature. Therefore, 




( 20 ) 


From substitution of equations (19) and (20) in equation (11) 
letting I = L, the effective length-diameter ratio of a flow 
passage can be expressed as 


'e/eff 


b.022\ 

0.023/ 






( 21 ) 
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In order to ascertain the validity of the assuiirption expressed 
hy eqixation (20)^ a con 5 )arison was made^ for several experimental 
runs, of the pressure drops measured in the investigation of refer- 
ence 6 with the pressure drops obtained from the charts using equa- 


passage. Inasmuch as in the investigation the wall temperature 
varied with length along the tube, the integrated average value 
of was used in the calculations* The results of the pressure - 
drop cQurparison as well as of a comparison of the cal ciliated and 
measured temperature -rise values are presented in table I. 

The good e^eement between the measured and the calculated 
results indicated in the table substantiates the genei^al validity 
of equations (19) and (20) and is a check on the accuracy of the 
charts and assuDiptlons involved therein. In addition, the results 
of the foregoing comparison show that, althou^ strictly applicable ^ 
for the case of constant wall teniperature, the charts can be used, 
within limits, to handle the case of nonuniform wall -temperature 
distribution along the passage length throu^ the use of the 
average wall ten 5 )erature. A measure of the nonunlformlty of wall- 
teniperature distribution is given by the 2 ?atio of the difference 
between max i m u m wall temperature and entrance -air temperature to 
the difference between average wall temperature and entrance -air 
temperature. For the runs used In the foregoing con 5 )arlson this 
ratio was approximately 1.2. 

At the relatively low Eeynolds number of 53, 600^ the agreement 
is not as close as at the higher Eeynolds numbers. This result is 
attributed to the laminar boundary- layer regime at the tube 
entiance which, because of the smooth well-rounded entry of the 
test tube, occupies an extensive portion of the tube length even 
for Eeynolds numbers greatly in excess of the critical value. The 
extent of the entrance leiminar regime is reduced with Increase in 
Eeynolds number. 


The method of determining f for a flow paesace is 

\®e/eff 

Illustrated later by means of an example. 

Flow in passeige with sudden increase in cross section. - In 
reference 5, the momentum equation describing the flow of a com- 
prassible fluid across a sudden enlargement of cross section is 
written in teims of dimensionless-flow parameters for convenient 
application to the case where the flow in the passage Just upstream 


tion (21) for evaluating the 



of the experimental flow 
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of the area enlargement is critical or supersonic. In heat- 
exchanger operation, however, the flow is generally subsonic; the 
general method of reference 5 is herein applied to this specific 
case. 

For subsonic flow across a sudden enlargement (from area A^ 
to area A 2 ), application of the principle of conservation of 
momentum at a section just downstream of the enlargement and 'at 
a section located a sufficient distance downstream of the enlarge- 
ment that uniform flow again exists gives: 


mVj^ + PpAg = ^2'^^ 


( 22 ) 


where 

subscript 1 denotes conditions at the small area A]^ 

subscript 2 denotes conditions at the large area A 2 

Because the mass flow and total temperature of the fluid are 
constant across the enlargement, equation (22) can be written 


/ mV+pA \ ^ mVi-)-PiA 2 ^ / mV+pA \ ^ 

m^gRTi 

or, more conveniently, when noting that f = 


Pi (A2 ~Ai) 
mVgRTi 

A2^ 


(23) 



The use of equation (24) for calculation of the compressible- 
flow changes across an enlargement for subsonic flow is illus- 
trated later by means of an example. 


DESCRIPTION OF CHARTS 

In figure 1, the ratio of the total momentum at any 

(mV+pA)^ 

station to the total momentum at the reference station ^mV+pAT 
is plotted against the reciprocal of the total -momentum 
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parameter at the reference station 

of effective length-diameter ratio 
_x 
J>. 


m VgRT 
mV+pA Jo 


for a range of values 


(De)eff* 


The ratio 


is taken for the distance between the reference station and 

’e/eff 

the station x in the flow system. Within each figure, the locus 
of all flow conditions in a given flow passage has one value 
of the abscissa. Hence, the points on a figure representing all 
flow conditions in a pass6ige lie on a common vertical line. Each 

curve of constant ^ figure 1 also represents the locus 


of points at which 
. T. ■ 

for 


(t". 


eff T 

— is a constant. Figure 1(a), which is 

■ X 


= 5.0, covers a range of 


a lent to a range of 


T 


'e/eff 


from 0 to 32.1 equlv- 


from 5.0 to 3.0; figure 1(b), which is for 


— ) =3.0, covers a range of 


'e/eff 


from 0 to 50.65 equivalent 


to a raj3ge of 


from 3.0 to 2.0; figure 1(c), which is for 


— ] = 2.0, covers a reinge of 

T/0 


T, 


■=— ) from 0 to 71.4 equivalent 
^e ^f f 


to a range of from 2.0 to 1.5; eind figure 1(d), which is for 


T, 


= 1.^ covers a range of 

T.. 


{%). 


from 0 to 208 equivalent to 
eff T„ 


a range of from 1.5 to 1.11. Thus a total range of — from 

5.0 to 1.11 is covered in figure 1. 

Inasmuch as the bottom curve in each of figures 1(a) 
to 1(c) is used in proceeding from one figure to another when 
required in solution of a problem, it is dashed and is identified 

with its appropriate value of . Ihese values are -^ = 3.0 
obtained at = 52.1 in figure l(a) , ^ = 2.0 obtained 




- (^’ 

fie-, 1(0). At (I-) 




, = 50.65 in figure 1(b), and -;^ = 1.5 obtained at 


= 0, -i^, of course. 


equals tlie value of 


(— ) 

\ T/q 


Because tlie value of corresponding to maximum 


for which the figure is presented. 


Jeft 
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on a figure is identical witii the reference value of 


(De)eff 


^ on the 


in figures 1(a) 


succeeding figure, the maximum value of 

to 1(c) represents the effective length -diameter ratio "between the 
reference station for which the figure is presented and the reference 
station in the succeeding figure. For example, from figure l("b). 


© 
(^1 = 


value between the reference station defined by 


h) = 


2.0 is 


3.0 and the reference station defined by 

'0 

50.65. The curve representing the choke limit for the flow is 
given in each figure. A table is presented in figures 1(a) to 1(c) 
that is useful, as later shown, in proceeding from one figure to 
another where roguired in the solution of a problem. 


Figures 2 and 3 relate for air (7 = 1.400) the total-pressure 


parameter 


PA 


and the static -pressure parameter 




respectively, to the total-momentum parameter ” The varla- 

m 

tlon of and T in a flow passage can therefore be trans- 

m 

lated Into the variation of static or total pressure in the flow 
passage. 


OEEEIATIONS INVDL7ED IN USE OF CHARTS 

The operations involved in the use of the charts for solution 
of a heat -exchanger flow problem are briefly outlined and illus- 
trated in the following exan^jles with the aid of figures 4 and 5, 
Figure 4 presents a schematic diagram of the flow passage of exam- 
ples I, H, and HI and indicates the passage (solid lines) as a 
segnent of a generalized flow system shown by the dashed lines. Ohe 
relative positions of the various stations in the generalized flow 
system that enter in iiie problem of example I are indicated. Bie 
detailed steps Involved in the use of the charts for solution of 
the problem of example I are traced in figures 5(a), 5(b), and 5(c). 
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Eiample I - Use of Chairts 

Consider the turbulent flow of air through a smooth tube 
wherein heat is added to the air at constetnt wall temperature. The 
following flow and heating conditions are assumed to exist; 


(1) Equivalent diameter of tube, Dq, (ft) 0.0833 

(2) Length of tube, L, (ft) 7 

(3) Mass flow of air throu^ tube, m, slugs/(sec) • * 0.00373 

(4) Total tenq)erature of air at tube entrance, ^en> 

^ 610 

(5) Total pressure of air at tube entrance, Pen^ 

(lb)/(sq ft) 2160 

(6) Q?ube-wall temperature, °E.. 2000 


Deteznlne: 

(a) Total pressiure and total temperature of air at tube exit 

(b) Static pressiire of air at tube entrance and exit 

The method of successive approximations is required to calcu- 
late the average value of ( for the tube passage that is 

\®e/eff 

consistent with the average fluid tempeiature in the passeige. 
Inasmuch as the principal purpose of this example is to illustrate 

the chart operations, the evaluation of al'bliou^ the 

first step in obtaining solution of the problem, is pi^sented later 
in exanq)le II. 


(7) The results of example II give 


'eff 


(i). 

(8) From items (4) and (6) 

2000 

V T / 


= 103.0 


en 


610 


= 3.279 




(' ^ 
(o| 


(9) From equation (15) and items (7) and (8) 


J 
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^5 ^ = 1 . (3.279 -1) (l - .-0.00568X103.0) . 

\(j\ 


2.010 


“■en 


(lO) From items (l) and (3) 




m 

pv = 7 = 


0.00373 


A It 


I X (0.0833)^ 


= 0.6837 


(11) From items (4), (5), and (lO) for E = 53.35 for air 

2160 

0.6837 V32.2 X 53.35 X 610 






= 3.087 


m 


(12) From item (11) and figvire 2 



3.252 


Deteiml, nation of initial chart. - Items (8) and (12) specify 
the entrance conditions of the tube required for entering the charts. 
The point on the charts corresponding to the entrance conditions of 
the tube is referred to as the starting point and the chart contain- 
ing the starting point as the initial chart. The initial chart is 

simply determined by the value of ( -~) as illustrated in the 

\ T /en 


following item . 


(13) From item (8), 
T 

covers the range ~ = 

T 



5.0 to 3 


= 3.279; hence, figure 1(a) -vrtilch 
.0 is the initial chart. 


Inasmuch as consideration of several reference stations in 
the flow system may be required in solution of a problem, the ref- 
erence station for which the initial cheirt is presented and all 
flow conditions at that reference station are designated by the 

subscript 0,1; for example , [ ■^) is equal to 5.0, the refer- 

\ T /o,i 

ence value on the initial chart. As shown in figure 4, reference 
station 0,1 in the flow system is located upstream of the flow- 
passage entrance (station en) . 
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Location of starting point on initial chart. - 


(14) 


en 

D. 


e/eff 


1b calculated fixan equation (14) for 


^)en " \^/0 


on the initial chart (designated 


[~^J equal to 5.0 from item (13)). 

^ ^ 0,1 


'^) =176.1 loge 

/eff 


5.0 - 1.0 


5.0 1- 


3.279 


= 24.76 


As shown in figure 4, \~^~] represents the effective 

\ ®e /eff 

length- diameter ratio between stations 0,1 and en, that is, the 

Ti 


w 


effective length-diameter ratio required to effect a change in 
from 5.0 to 3.279. 


(15) The ordinate-to-abscissa ratio on the initial chart for 


X = X, 


en 


is 


(mV+pA)en / / mV^ \ 
(mV+pA)^^^ / V“V+pA Jq^ 2. 


which, because is constant, can be rewritten as 


/ mV+pA \ 

VmV^Ven V ^^^4n ^ "^0,1 


This ratio^which represents the slope of a straight line through 
the origin of coordinates in the initial chart, is evaluated from 
items (8), (12), and (13) as 


3.252 


i 


5.0 

3.279 


4.013 
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Intersection of the straight line with slope 4.013 and the 

interpolated curve for = 24.76 (item (14)) 

\I>e/eff \®e/eff 

locates the starting point on the Initial chart. Inasmuch as the 
origin of coordinates is not on the charts, the required slant line 
is obtained hy drawing a strai^t line through the points 
(ordinate = 1.00, ahsclssa = 1/4.013) and (ordinate = 0.94, 
abscissa = 0.94/4.013). The foregoing construction procediire for 
location of the starting point is illustrated in figure 5(a) where 
the starting point is designated as point A. 

(16) The ordinate value for the starting point 

^Ich is read from figure 1(a) as 0.990. 


(17) The abscissa value for the starting point 
which is read from figure 1(a) as 0.2467. 

Determination of terminal chart. - The terminal chart is that 
chart containing the point corresponding to the exit conditions of 
the flow passage. This point is herein referred to as the end 
point and is designated as point D in figure 5(c). Hie te rm i n al 
chart is determined as Illustrated in items (18) and (19). 

(18) The effective length-diameter ratio between reference 
station 0^1 and the end of the tube passage (station ei in fig. 4) 

is — ) , which is simply the sum of ( -—] (item (14)) 

V ^e Jeff \ ^e Jeff 

and the | of the tube passage (item (7)). 

V^e/eff 


is 


(mV+pA)en 

(mV+pA)Q^- 


is 


m 


VgRT 


mV+pA 


0,1 


. /eff 


24.76 + 103.0 = 127.76 


(19) The table presented on the initial chart (fig. 1(a), in 
this case) specifies the terminal chart as fig. 1(c) for the value 


^en"*"^ 


'e /eff 


equal to 127.76. 


Transfer from initial to terminal chart. - Inasmuch as the 
terminal chart is not the same as the initial chart, solution of 
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the problem requires transferring from chart to chart until the 
terminal chart is ^tered. In this case, it is necessary to trans- 
fer successively from figure 1(a) to figure 1(b) and finally to fig- 
ure 1(c). In the performance of this transfer, the abscissa values 
required to enter fig;iree 1(b) emd 1(c) are successively determined 
as follows; 

(20) On the initial chart, the ordinate corresponding to 

) along the curve for marimvnn ( ) is equal to 

0,1 \®e/eff 



(mV+pA)n 2 

I ' . \ * where subscript 0,2 denotes the reference station at which 

(mV+pA)o^l 


_w 

T 


V 

T 




= 1—1 = the value of -ts- at maximum ( ^ 


0.2 


on the 


e/eff 


Initial chart. From figure 1(a) for item (17) and for 
= 32.1 (maximum value). 


e/eff 


(mV-HpA)o 2 _ 

(mV+pA)Q^ 


= 0.9862 


•<diere 


= 3.0 

/ 0,2 

As indicated in figure (4), station 0,2 is the next successive ref- 
erence station downstream of station 0,1. The point on the charts 
cori'esponding to the foregoing conditions is designated as point B 
in figure 5(a) . 

(21) The reciprocal of the total-momentum parameter at refer- 
ence station 0,2 is calculated as 


m VgRr \ _ f m_V^®'\ (°»V^-pA) 


mV+pA Jq 2 \mV+pA 


0,1 
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•which from items (17) and (20) for 


I2f 

T 


3 

T 


0,2 


= 3.0 (item (20)), gi-ros 


0.2467 

0.9862 




0,1 


3229 


= 5.0 (item (13)) and 


The chart for =(~^] =3.0 (fig. iCb)) is entered 

\T/q \T/o,2 

with an ahscissa value equal to 0.3229 as indicated in figure 5(h). 
Inasmuch as figure 1(h) is not the teminal chart, "bhe procedure of 
items (20) and (21) must he repeated. 


ponding to 


(22) On the chart for 
( m 


T, 


the ordinate corre- 


T/O " VT/0,2 

along the curve for maximum 


is 


eff 


\^mV+pA 70^2 
(mV+pA)o 3 

equal to -r _ where subscript 0,3 denotes the reference station 


(mV+pA)Q'2 


at which I 


the value of at maximum 


(®e)eff 


on the chart. From figure 1(h) for item (21) and for 

/-£\ = 50.65, 

VDe^ff 


(mV+pA)o 3 

(n>v+p^)o;2 ^ 


where 






As indicated in figinre 4, station 0,3 is the next successive 
reference station downstream of station 0,2. The point on the 
charts corresponding to the foregoing conditions is designated as 
point C in figure 5(h) • 
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(23) The reciprocal of the total-momentum parameter at refer- 
ence station 0,3 is calculated as 


W^\ //^\ 7t\ 

/0,3 ■ Jo,Z ("V+pAJo’j yV I^,2 


Which, from items (21) and (22) for =3.0 (item (20)) and 

0^2 


~^) = 2.0 (item (22)), gives 


0,3 


0.3229 

0.9621 


_ n 

V 2.0 “ 


4111 


The chart for 


3 

T 


= y-^j = 2.0 (fig. 1(c)) is entered with an 
^ '0,Z 


abscissa value equal to 0.4111 as indicated in figure 5(c). 
Inasmuch as figure 1(c) is the terminal chart (item (19)), the 
transfer procedure is completed. 

Location of end point on terminal chart. - 

(24) The table presented on the initial chart (fig. 1(a), in 

this case) specifies Z = 82.75 for (— — ) = 127.76 

(item (18)). 


De Jeff 


(25) The end point is located on the curve in the terminal 
chart for 


(l>e)eff 




'e /eff 


which, from items (18) and (25), gives 


on terminal chart = 127.76 - 82.75 = 45.01 


'e/eff 


(26) Items (23)and (25) fix the location of the end i>olnt on 
the terminal chart (fig. 1(c)); item (23) specifies the abscissa 
value for the end point as 0.4111 and item (25) specifies the 
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( =r- 1 value for the end point as 45.01. This procedure is illus 
\I^e/eff 

trated in figure 5(c) where the end point is designated as point D. 


Evaluation of 


/ mY+pA \ 

\idV 


(27) Note that 3 charts, and thus 3 reference stations are 
involved in this problem for location of the end point. In this 
case, the ordinate value for the end point is 


(mV + pA)q^ 

(»V+pA)o ^3 


which, from item (26), 

(mV+pA)^^; 

^28) (mY+pA)Q^3^ 


is read as 0.9463 in figure 1(c). 
is calculated as 


(mV+pA)^^ (mV+pA)Q^2 (mY+pA)^^^ (mV+pA)^^ 

(mV+pA)^^^ (mV+pA)^^^ (mV+pA)Q^2 (inV+pA)^^^ 


^diich, from items (20), (22), and (27), is 


0.9862 X 0.9621 X 0.9463 = 0.8979 


(29) 


(mV+pA)^^. 

(■»V+pA),^ 


is calculated as 


(mV+pA)^^. (mV+pA)^^ (mV+pA)o^i 

(mV+pA)^^ " (“V+pA)q^j^ (mV+pA)^^ 


which, from items (16) and (28), gives 


I 


(mV+PA)ei ^ 0.8979 

(mV+pA) 0.990 
^ en 


0.9069 


1004 


NACA EM No. E8G23 


27 


(30) The total-momentian parameter at the exit of the flow 
passage Is then 


/ mV+pA \ / mV+pA \ , __ 

1/ ^ei 

6H 


en 


which, from items (9), (12), and (29), gives 

3,252 X 0.9069 


/ m7+ pA, 

V m VS^^ex 


V 2 .OIO 


= 2.080 


Final results. 


(31) From item (30) and figure 2(a) 


PA 


m 




= 1.778 


ex 


(32) From items (9), (11), and (31) 


•ex 

^eu 


m 




010 = 0.817 


ex " ' en 

( 33 ) From item (12) and figure 3(h) 

£A 


m 




= 2.913 


en 


( 34 ) From item (30) and figure 3 


ym 


= 1.420 


®V^;ex 

( 35 ) From items (5) and (32) 

= 0.817 X 2160 = 1765 (lh)/tsq[ ft) 

(36) From items (11) and (33) 


en 


3.087 

2.913 


= 1.059 
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so that, from item (5) 

I’ea = = 2040 (Ibj/faq ft) 


(37) From items (31) and (34) 

/P\ 1.778 

\vL = 1.420 


1.252 


so that, from item (35) 

“ 1.252 ■ (l^y(s^ 

(38) From items (4) and (9) 

Tqj. = 2.010 X 610 = 1226° E 


Example II - Evaluation of (fr* ) 

The method of successive approximation must he used in order to 

obtain consistent values of { and average fluid temperatiire. 

\i'e/eff 

The variations of specific heat Op and absolute viscosity (x of 
air with temperattire, as obtained from reference 8, are presented in 
figure 6 to aid in this deteimination. 

(39) From figure 6 for T^ = 2000° E ^item (6)) 

%,v = 0.2775 

= 30.3 X 10"® 

T^y = 750° E as a first approximation in which case 


°p,av 


(40) Take 
from figure 6 


= 0.2425 


1004 


^001 . 


(41) From eq.mtion (21) and items (l), (2), (6), (10), (39), and (40), noting that 
equals 12.3 X lO"® (SYMBOIS) 

4L ■ (fei) (sr (as)(s:s5)"'‘ fe) 


=89.1 first approximation 
(42) From equation (15) and items (8) and (41) 


en 


= 1 + (3.279- 1) = 1.9052 first approximation 


(43) From items (4) and (42) 

T. 


1 + 


^av “ ^en 


'•ex > 

'^en 


= 610 


^1 +1.9052^ 


886° E second approximation 


in which case from figure 6, Cp^av = 0.2455 

(44) From equation (21) and items (l), (2), (6), (10), (39), and (43) 

4L ■ (SI) (S8) (“S)“ kkd 


M2. 3X10 

= 100.6 second approximation 
(45) From equation (15) and items (8) and (44) 


CO 
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^ejc 1 ^ P7Q ^-0. 00568X100. 6\ ^ ^ 

= — = 1 + (3.279- l; e )= 1.993 second approximation 

■‘•en 


(46) From items (4) eind (45) 

Ta^ = 610 

in viaich case, = 0.2463 


= 913° R third approximation 


(47) From equation (21) and items (l), (2), (6), (10), (39), and (46) 

® (^ 33 - 

= 102.7 third approximation 

(48) From equation (15) and items (8) and (47) 


*1 , 07 Q T \ —0.00568X102.*^ — r\r\'7 4-i.4«;a • 4.4 

7 ;; = 1 + (3.279- 1; \1 - e y = 2.007 third approximation 


T. 


en 


(49) Repetition of the foregoing procedure gives for the fourth approximation 

Tav = 917° R 


r^) = 103.0 

I>e/eff 


•■ex 


= 2.010 


■•en 


o 
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wbich are practically identical with the third approximation values. 

T 

The values of Tg^y, , and 7=;^ given in item (49) 

\"e/eff ■‘■en 

represent the convergent values. It is evident that the third 
approximation values could have been taken as the convergent values 
with negligible error; in general^ the third approximation values 
are sufficiently accurate. 


Example HI - Determination of Exit End Losses 

Assume that the flow at the exit of the tube of exan5)le I 
suddenly discharged into a duct having twice the cross-sectional 
area of the tube (in which case, f = 0.50), 

Detenaine : ^ 

(a) Static and total pressures at the duct cross section where 
imifom flow is reestablished (referred to as station d in fig. 4) . 

(b) Drop in static and total pressures from entrance of tube 
of example I to duct cross section where uniform flow is reestab- 
lished (in fig. 4, fivDm station en to station d). 

(1) From the results of exaarole I (items (30), (31), and (34)), 
the flow parameters at the tube exit are: 


.2.080 



1.778 



1.420 


(2) From equation (24) and item (l) of this example for f = 0.50, 
noting tlrnt in figure 4 the small area is at station ex and the 
large area at station d 
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m7+^\ ^ = 3.500 

ttia/^/, W.50 y 


V^ya 

( 3 ) From item (2) and figures 2 and 3 


PA 


m VsKT 


= 3.348 


= 3.189 

m V^/d 


’•ei 


( 4 ) From items (l) and (3) noting that = f = 0.50 emd 

that m and T are constant across the enlargement 


PA 


62: \m VsKP/d 


' m 


^ X 0.50 = 0.942 


PA 


ex 


Ad ■ 1.778 


- / pA \ [mV^^j _ei ^ 3J^ q^ 5 q ^ 1.123 

^ex \inV^j^ ^ ^ /ex ^d 

( 5 ) From example I, P^^ = 1765 pounds per sqxiare foot and 
Pg^ = 1410 pounds per square foot, which from item (4) gives 

P^ = 0.942 X 1765 = 1663 (lh)/(sq ft) 

p^ = 1.123 X 1410 = 1583 (lb)/(sq ft) 


(6) The drop in total pressure from station en to station d 
in figure 4 is, from item (5) of example I and from item (5) of 
this example, 

dPen_d = 2160 - 1663 = 497 (lh)/(sq ft) 

If the flow process at the tube entrance were isentropic, the value 
of total-pressure drop given by item (6) would also be obtained 
across the entire tube, including all end losses. 
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(7) The drop in static pressure from station en to station d 
in figure 4 is, fran item (36) of example I eind from item (5) of 
this example, 


APen_d = 2040 - 1583 = 457 (lh)/(s(i ft) 


The static-pressure drop obtained at the tube entrance for an 
isentropic flow process can be calculated by standard equations and 
added to the value of pressure drop given by item (7) to obtain the 
over-all static-pressure drop across the tube, including all end 
losses* 
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TABLE I - COMPARISON OF CALCULATED AND MEASURED^ RESULTS 


Average 

Tv 

(°E) 

Entrance conditions 

Exit 

M 

Msasured. results 

Calculated results 

Percentage difference 

Reynolds 

number 

P 

(It/sq ft) 

T 

(^) 

M 

Pen 

Tex 

T 

en 

Ap 

(ib/sq ft) 

AT 

(°E) 

Pen 

T 

ex 

T 

en 

Ap 

(It/sq ft) 

AT 

(°R) 

in 

Ap 

(Ib/sq ft) 

AT 

(°E) 

1672 

184,000 

4527 

528 

0.37 

0.84 

0.551 

1.710 

2031 

375 

0.544 

1.699 

2065 

369 

1.7 

-1.6 

1072 

155,000 

3682 

554 

.41 

.70 

.639 

1.359 

1329 

199 

.657 

1.348 

1264 

193 

-4.9 

-3.0 

861 

219,000 

4458 

540 

.46 

.87 

.568 

1.218 

1926 

118 

.553 

1.218 

1993 

118 

3.5 

0 

1670 

263,000 

6186 

530 

.39 

1.00 

.464 

1.651 

3315 

345 

.458 

1.645 

3353 

342 

■^1.2 

-1.0 






choke 











1752 

53,600 

2417 

527 

.21 

.32 

.881 

1.967 

288 

510 

.889 

1.939 

267 

495 

-7.1 

-3.0 


^■Data from reference 6. 

^Calculated choking length less than actual choking length hy 2.6 percent. 



04 

VJl 


NACA RM No. E8G23 


Ratio or total momentum 


c^ 





.16 .20 .24 .28 .32 .36 .40 .44 .48 .52 .56 


Reciprocal of total-momentum parameter at reference station. 


(a) 


Reference station characterized by 



5.0. 



Figure 1. - Variation of total momentum with distance 




eff 


at the reference station. Tiirbulent flow, air (y = 1.400), 
(A 7 - by 22-inch print of this figure is attached.) 


as a function of the total-momentum parameter 
heat addition at constant wall temperature. 
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Figure 1. - Continued. Variation of total momentum with distance ( ^) as a function of the total-momentum 

j parameter at the reference station. Turbulent flow, air (y^ = 1.400), heat addition at constant wall 

temperature. (A 9- by 22-lnch print of this figure is attached.) 
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Ratio of total momentuin, 1 
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Reciprocal of total-momentiam parameter at reference station. 


^mV+pA 
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Figure 1. - Continued* Variation of total momentum with distance ( ^ ) as a function of the total -momentum 

\ ^/eff 

parameter at the reference station. Turbulent flow, air (y - 1.400), heat addition at constant wall 
temperature. (A 10- by 22 -inch print of this figure is attached.) 
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Figure 1. - Concluded* Variation of total momentum with distance as a function of the total- 

Uft . 

momentum parameter at the reference station* Turbulent flow, air (y » 1.400), heat addition at constant 
wall temperature* (A ib- by 22-inch print of this figure is attached*) 
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(a) Range of — 1.852 to 2,42. 

m VgRT 

Figure 2, - Relation between total-pressure parameter and total-ncmi-ntura parameter for air 

(Y = 1*400). 
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Flgiare 2. - Continued, Relation between total-prese\n»e parameter and total -momentum para- 
meter for air (y ■ 1«400)« 


mVS^ 


42 


NACA RM No. E8G23 



mV ♦ pA 

mVgRT” 

(c) Range of 3,75 to 5.15. 
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Figure 2. - Continued. Relation between total-pressure parameter and total -momentum para- 
meter for air (y =» 1.400). 
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Flgizre 2« - Concluded. Relation between total -pressure parameter and total^omentiim para< 

meter for air (y =* 1,400). 
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(a) Range of — 1.852 to 2.32. 

IB 

• “ Relation between atatlc-pressure parameter and total^^oonientum parameter for air (y » 1.400) • 
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(b) Range of 2.30 to 3.85. 

m VgRT 

Pigitre 3. - Continued. Relation between statlc-preseiire parameter and total -momentian parameter for air 

(y « 1.400). 
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Figure 3* - Continued* Relation between e tatlo-preeiiure parameter and total-momentum parameter for air 

(y « 1.400). 
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• Relation between scatlc-pressiire parameter and total -momentum parameter for air 
(Y = 1.400). 
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Figure 4# - Schematic diagram of flow passage of examples I, II, and III showing relative positions of 

various stations in given flow system. 


E8G23 



(a) Figure 1(a), (^) = 


Figiire 5* - Use of figure 1 in solution of problem of example I. 
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(b) Plgiore 1(b), (“»f) “ 3.O. 
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Figure 5, - Continued. Use of figure 1 in solution of problem of example I. 
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Figure 6. - Variation of absolute viscosity [i and specific heat Cp of 

air with temperature. 


